One dimensional scattering of a two body interacting 

system by an infinite wall 



A.M. Moro, J. A. Caballero and J. Gomez-Camacho 

Departamento de FAMN. Facultad de Fisica. Universidad de Sevilla. 
Apdo. 1065, 41080 Sevilla (Spain) 



Abstract 

The one-dimensional scattering of a two body interacting system by an infinite wall is studied in 
a quantum-mechanical framework. This problem contains some of the dynamical features present 
in the collision of atomic, molecular and nuclear systems. The scattering problem is solved exactly, 
for the case of a harmonic interaction between the fragments. The exact result is used to assess the 
validity of two different approximations to the scattering process. The adiabatic approximation, 
which considers that the relative co-ordinate is frozen during the scattering process, is found to 
be inadequate for this problem. The uncorrelated scattering approximation, which neglects the 
correlation between the fragments, gives results in accordance with the exact calculations when the 
scattering energy is high compared to the oscillator parameter. 

PACS numbers: 24.10.Eq;03.65.Nk;3.80.+r;24.50.+g 
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I. INTRODUCTION 



Recently, important efforts in the fields of molecular, atomic and nuclear physics have 
been devoted to the analysis of collision processes involving composite quantum systems. 
Despite the peculiarities of the different fields, the theoretical description of the collision 
has several common features. R-matrix theory [l] is used for re-arrangement collisions, 
when several electrons or nucleons may be exchanged between the colliding systems. Cou- 
pled channels calculations jl, 2| are used to describe excitation and dissociation in atomic, 
molecular and nuclear systems. The "adiabatic" or "sudden" approximation is often invoked 
in molecular and nuclear collisions, because it simplifies significantly the description of the 
scattering process, by considering that some of the relevant co-ordinates are effectively frozen 
during the scattering process. Also, the description of the collision of atomic and molecular 
beams with surfaces requires coupled channels descriptions, but the combined difficulties of 
the atom-surface interactions and molecular vibrations makes the "sudden" approximation 
almost essential for its solution |2|. 

In the case of atomic and molecular physics, one-dimensional atom-molecule collision reac- 
tions have been studied in detail by solving the Schrodinger equation In particular, colli- 
sions between an atom and a diatomic molecule represented by harmonic [4| , an-harmonic Q 
and Morse oscillators j^l have been analysed with different degrees of approximation. Diverse 
computational methods have been also implemented to study three-dimensional molecular 
and atomic collisions Moreover, in some recent papers algebraic approaches have been 
proposed for describing one- and three-dimensional atom- molecule collision processes js, 9|. 

In the field of nuclear physics, much of the interest in recent years has been focused on the 
study of the properties of halo nuclei, weakly bound systems characterised by the existence 
of one or two particles (generally neutrons) with a high probability of being at distances 



larger than the typical nuclear radius. Different approaches have 
of reactions involving halo nuclei. The adiabatic approximation 



aeen used in the analysis 
10-13] assumes that, for 



sufficiently high scattering energies, the internal Hamiltonian is accurately represented by 



its corresponding eigenvalue for the ground state. The sudden approximation |14H17] re- 
lies on two main assumptions: i) the impulse approximation, i.e., the multiple scattering 
expansion for the T-matrix is approximated by the sum on the individual T-matrices for 
the scattering of the separated constituents, and ii) only one of the particles of the pro- 
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jectile interacts with the target. In the high energy regime, the Glauber approach 18], 
which combine eikonal dynamics with the adiabatic approximation, provides a simple tool 
to analyse reactions involving halo nuclei. We have recently developed an alternative ap- 
proach to the description of weakl y b ound systems. The approach, called "Uncorrelated 



Scattering Approximation" (USA) 19|, |20j, is based on the fact that, for a weakly bound 



projectile, the correlations between the constituents are weak and so, to some extent, they 
are expected to evolve independently in the strong field of a heavy target. Thus, the three 
body S-matrix can be expressed in terms of the individual two body S-matrices for the 
scattering of the constituents. In this framework, the scattering observables of the process 
are given in terms of two body constituent-target observables. This model has been applied 
to describe elastic scattering and break-up of deuterons on heavy targets, with encouraging 
results. The Uncorrelated Scattering Approximation has certain relation with the R-matrix 
approach. In a R-matrix calculation of deuteron scattering, the wave-function within the 
range of the target interaction is given in terms of products of single-particle (protons and 
neutrons) wave-functions with the proper boundary conditions, provided the interaction 
between the fragments of the projectile is neglected. These single-particle wave-functions 
should be matched with the proper asymptotic wave-functions. In the USA calculation, the 
incident wave is expanded in terms of products of fragment-target wave-functions, which 
then scatter independently from the target. 

The objective of this work is to investigate the scattering of a composite system from a 
target with interactions which have a very short range compared not only with the size of 
the system, but also with the associated wave-length of the projectiles. For that purpose, we 
consider the case of two particles, interacting through a harmonic oscillator potential, which 
collide with an infinite wall. We develop two different methods to solve the problem exactly, 
obtaining the wave-function as well as the S- matrix, or reflection coefficients, which give the 
probability amplitudes for the excitation of the different oscillator states. We compare the 
exact results with the ones obtained using the adiabatic approximation, and the uncorrelated 
scattering approximation. 

The model discussed here does not pretend to be a realistic representation of any specific 
molecular, atomic or nuclear system. However, it has the virtue that the only length scale is 
the oscillator length ao, while the only energy scale is the oscillator parameter Tiu. Then, the 
results obtained, which are expressed in terms of rja^ and E/hu), can be applied in principle 
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FIG. 1: Schematic plot representing the scattering problem in the (r, R) plane. The waves represent 
the SCLS associated with a basis with N = 8 HO states. Each one of these SCLS has been plotted 
at the center of mass distance at which it is supposed to be scattered by the wall, according to the 
formalism described in section Hi Bl 

to arbitrary energy and length scales, which may be nuclear, atomic or molecular. This fact 
makes the model attractive as a bench-mark to test the validity of the approximations which 
are used for the description of composite systems. 

One aspect of the model which seems odd is the infinite nature of the harmonic oscillator 
interaction between the fragments. This would prevent dissociation in the case of molecules, 
ionization in the case of atoms or break-up for nuclei. On the other hand, the harmonic 
oscillator basis is complete, and thus the results of the model incorporate effects on the 
scattering due to coupling to all possible, open or closed, states. 

The paper is organised as follows: in section II we describe the problem to be solved. Two 
different approaches to extract the exact solutions are developed. In section III we derive 
the S-matrix in the adiabatic model approach. In section IV we present the Uncorrelated 
Scattering Approximation. In section V we apply it to evaluate the scattering coefficients 
and compare with the exact solution. Section VI is devoted to summary and conclusions. 

II. STATEMENT OF THE PROBLEM AND EXACT SOLUTION 

Let us consider the scattering of a two-particle bound system by the short range potential 
due to a heavy target, placed at the origin of coordinates. 

The Hamiltonian for the system may be written in terms of the coordinates of the two 
particles (r 1; r 2 ) and their momenta (pi, P2) as 
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FIG. 2: S-matrix elements versus the scattering energy for the HO states n =0, 2 and 4. The solid 
and dashed lines correspond to the exact calculation using the methods of sections III Al and III Bl 
respectively. A basis space with 15 even parity HO states was used for the calculations. 



H = h x + h 2 + v(n - r 2 ) (1) 

fc. = |: + ".M (2) 

where are the masses of the constituents, Vifa) is the potential exerted by the target on 
each particle and v(ri — r%) represents the binding interaction. 

Alternatively, it can be expressed in terms of the relative (r) and centre of mass coordinate 
(R): 



H = H r + H R + Vi(R + a 2 r) + v 2 (R - a x r) 
A = f +t,(r) 



(3) 
(4) 
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where M is the total mass, aj = rrii/M (i = 1, 2), and H r is the internal Hamiltonian. Here, 
[i represents the reduced mass of the two-particle system and r = r\ — r 2 is the relative 
coordinate. 

We refer to the specific case in which fj(rj) corresponds to an infinite wall potential, i.e.: 

, n > 

vM) = { (6) 

OO Ti < 

For negative values of r\ and r 2 the total wave function must vanish. For positive val- 
ues, the Hamiltonian comprises two terms: one associated with the internal motion of the 
projectile (H r ) and the other describing the centre of mass motion (Hr). 

We take v(r) to be a harmonic interaction. Therefore, the eigenfunctions for the internal 
Hamiltonian are given by 

0n(r)=A/^expf-^W(^) ; n = 0,l,... (7) 



where do = y %/ fioo is the oscillator length, M n a normalisation constant and H n the Hermite 
polynomial of order n. 

Denoting by iV the number of open channels, the total wave function will be then ex- 
panded in terms of eigenstates of the Hamiltonian as follows: 

i N—l o oo 

¥(r, R) = -^4>o(r)e- tKoR - £ ^4> m (r)e^ R - £ F Qm <p m (r)e~\ K -\ R . (8) 

V V m=0 V V m m= N 



In this expression K rn represents the centre of mass momentum associated with the in- 
ternal state m and v m = hK m /M its velocity. Energy conservation applied to the whole 
system leads to the constrain 

^ + , m = £, (9) 

where e m is the eigenvalue corresponding to the internal state <f) m (r). 

The first term in (jSJ) represents an incoming wave, normalised to unit flux, coming from 
R = +oo. The second term contains the set of scattered waves travelling in the positive 
R direction corresponding to open channels, i.e., e m < E. The coefficients S nm are the 
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FIG. 3: Projection of the total wave functions on the internal eigenstates n—0, 2, 4 and 6, for the 
scattering energy E = 5 %U). The solid line corresponds to the exact calculation by direct application 
of the boundary conditions, whereas the dashed line corresponds to the exact calculation in the 
CLS approach. 

S-matrix elements or, strictly speaking, reflection coefficients. Then, So m represents the 
amplitude probability of populating the state m during the collision starting from a wave 
function in its ground state. The last term in (jBj) contains the contribution to the wave 
function due to the (infinite) set of closed channels, for which e m > E. For these states, the 
associated momentum K m is a pure imaginary quantity, and the centre of mass motion is 
described by an exponential decaying behaviour. Therefore, they do not contribute to the 
asymptotic wave function and so they do not give direct contribution to the outgoing flux. 
Although they are usually ignored in practical calculations, the peculiarity of the infinite 
potential requires however the inclusion of these states in order to describe correctly the 
wave function in all the space. 

In this problem there is no transmitted wave due to the presence of the wall and so the 
S-matrix should fulfil unitarity. Therefore, the ingoing and outgoing flux must be equal: 
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N-l 

E \ S nm\ 2 = 1, (10) 
m=0 

where the sum extends to the set of open channels. 

The coefficients S nm and F nm are determined by imposing the boundary conditions at 
Ti = 0. For the infinite wall the wave function must vanish at r\ — and r 2 = which, in 
terms of r and R means 



\&(r,R = +a t r) =0 (r > 0), 

V(r,R = -a 2 r) =0 (r < 0). (11) 
In what follows we present two alternative methods to solve exactly this problem. 

A. Exact solution by direct application of the boundary conditions 

As described above the exact solution of the scattering problem of a two-body projectile 
by a rigid wall is accomplished by applying the boundary conditions (TTTj) to the general 
solution (jSJ). For simplicity we assume that the system is initially in its ground state. The 
first condition in (TTTI) leads to the following equation for the scattering coefficients: 

i N—l c oo 

lKoair Mr) ~ E ^e iKnair <t>n{r) - E F^ K ^ r (f>n{r) = ; r > 0. (12) 



e ^ 

V n =0 V Vn n=N 

The second condition in (TTTj) leads to a similar equation (with a>2 instead of ai) which holds 
for r < 0. In particular, we have performed calculations for the particular case of equal 
masses, i.e., a% = oli = ^, for which both equations are identical, due to the symmetry of 
the problem under the exchange of r\ and In the remaining of this section we restrict to 
this particular case. 

In order to transform eq. (Tl2|) into an ordinary algebraic equation, we require that 

1 2 

ty(r,R=-r) dr = 0, (13) 
which gives for the scattering coefficients the relation: 

J oo oo 

p _ E (Con°On + C()n a 0n) + E ^OnCl m b nm = 0, (14) 
^ n n,m 
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where 



c, 



On 



^S 0n ; n<N 



On 



n> N 



(15) 



and 



aon = / <P* (r)Mr)e liKn+Ko) i r ; n < N, (16) 

JO 

bum = / 0;(r)0 m (r)e^-^^ ; n , m < iV. (17) 

J 

The expressions for a Qn and 6 nm for n,m > N are obtained substituting ±iK n for — |i^ n |. 
Differentiating eq. (1141) with respect to Cq u we get the following linear system in the 
variables C 0n : 

oo 

"^2 Combnm = a* Qn , 71 = 0, . . . , OO. (18) 
m=0 

For practical calculations the resolution of this system requires the truncation of the sum 
at some finite value of m. It should be noticed that in order to achieve convergence for 
the S-matrix elements one is forced to include in the calculation several closed channels. 
Otherwise, the boundary conditions (fTTj) are not accurately fulfilled. 

B. Exact solution using a discrete basis 

In this section we present an alternative method to obtain the exact solution of the 
problem stated above. The method relies on the introduction of a new basis of states which 
are particular linear combinations of the internal wave functions. The new functions have 
the property of being highly localised in configuration space. As we shall see this peculiarity 
allows to apply more easily the boundary conditions. 

We start with a truncated basis of iV eigenstates for the internal Hamiltonian that we 
denote by {\Nn);n = 1, . . . , N — 1}. Thus, according to our previous notation, (r\Nn) = 
<Pn{r). 

In the appendix we show how these states can be decomposed in terms of configuration 
localised states (CLS) as 
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N 



\Nn) = J2 ( CLS ^ Ns\Nn) \CLS; Ns) (19) 



s=l 



In this expression the ket \CLS; Ns) represents a configuration localised state. Explicit 
expressions for the CLS associated with the HO wave functions can be found in the appendix 



and in Ref. (22). The function (r\CLS; Ns) has the property of being highly localised around 
r = r s , the s-th zero of the eigenfunction (r\N N). 

The problem is significantly simplified in the case of fragments of equal masses. Par- 
ity conservation guarantees that only states with the same parity as the incident wave 
function will be populated in the process. This allows to work in the subspace: {0„(r); 
n = 0, 2, . . . , N — 2} (for simplicity of the notation, and without loss of generality, we take 
N even). These functions are symmetric with respect to their natural variable, r. From this 
set of y states it is possible to construct, by means of a transformation similar to f TT9"j) . a 
set of symmetric configuration localised states (SCLS), which are also even functions with 
respect to the variable r. We denote this new set of states by {\SCLS; Ns); s = 1, . . . , y}, 
where the index s runs over the positive zeros of Hn(x). The details of its derivation are pre- 
sented in the appendix. The state (r\SCLS; Ns) has the property of being localised around 
the points r = r s = ±x s ao, where x s is the s— th positive zero of the Hermite polynomial 
T-Ln{x). In the treatment that follows we make extensive use of this remarkable signature. 

The boundary condition due to the wall requires that the total wave function vanishes 
for 7"i = and r 2 = or, in terms of coordinates r and R, along the lines R = r/2 (r > 0) 
and R = —r/2 (r < 0) in the (R, r) plane. At each value of R a reflected wave is generated, 
interfering with the other outgoing waves to construct the total scattered wave. This picture 
is simplified working in terms of SCLS. The part of the wave function associated to the state 
(r\SCLS; Ns) is peaked around r = ±r s and, therefore, it will be mainly scattered around 
R — R s — \r s \/2. In the limit case iV — > oo, (r\SCLS; Ns) becomes a delta function in 
r and the associated wave is exactly scattered at R — R s . Moreover, continuity of the 
wave function implies that a reflected wave, affected by a phase factor — exp(— 2iK R s ), is 
generated at this point. Obviously, this is not exactly our situation as, in practise, iV is 
always finite and so our localised states have a finite dispersion around r s . However, we can 
make this dispersion as small as required by increasing the number of states. 

Consequently, the process is considered as a distribution of localised states that are 
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reflected at some definite barriers in the R direction. The total scattered wave is given by 
the superposition of the reflected waves. To make the treatment clearer, we divide the R axis 
in — regions, delimited by the values of R s . We order the zeros of the Hermite polynomial 
%n in such way that R\ > R2 > . . . > Rn/2- Let us introduce an index i to label each region 
(i = 1, . . . , y), such that i = 1 corresponds to the asymptotic region, i.e., R > R\, before 
any localised state has been reflected. In this region the basis space associated with the 
internal motion is described in terms of the first y HO eigenfunctions with positive parity. 
Alternatively, it can be described in terms of y symmetric localised states. 

At R = R\ the SCLS corresponding to s = 1 is reflected and removed from the incident 
wave function, while the rest of SCLS remain unaltered. Therefore, in the region i = 2 
the basis is limited to the subspace spanned by the remaining y — 1 states (s = 2, y). 
Subsequently, our original set of states are no longer eigenstates of the Hamiltonian in this 
region. Instead, a new family of y — 1 eigenstates must be calculated, by diagonalizing the 
Hamiltonian in the basis constituted by the remaining y — 1 localised states. 

The method is schematically illustrated in Fig. [U The picture represents the scattering 
problem in the (r, R) plane. The total wave function travels in the R direction (vertical axis) 
and must vanish along the lines r\ = and r 2 = 0, which have been also plotted for reference. 
The case with N = 8 has been considered, in which the incoming wave is decomposed in 
a set of four SCLS. Each one of these SCLS is consider to scatter at a definite barrier in 
the R axis, labelled from Ri to R4. Note that they partially extend to the forbidden region 
(ri < 0), which is a consequence of the truncation of the HO basis. 

Let us specify explicitly the boundary condition at R = Ri. This barrier separates 
regions i and i + 1. The eigenstates corresponding to region i will be denoted by {|iV"W m); 
m — 0, iVW — 1} where iVW = y — i + 1 is the number of states in region i. The total 
wave function is then expanded in each region in terms of the corresponding eigenstates. 
For region i, corresponding to the interval Ri-\ > R > Ri, we have: 

\*®(E))= Yl (A^e- lK - R -B^e +lK - R )\N^m), (20) 

where A® (B$) are the coefficients of the incoming (outgoing) waves. Energy conservation 
requires, in analogy with Qj, 



11 



1 - 



eg 



o 
o 



CO 



J_ 0.1 - 



0.01 



~ i 1 1 1 1 1 1 1 r 



_i i u L 







Exact 

— Adiabatic 



5 10 
E (/zco units) 



15 



FIG. 4: Representation of the quantity 1 — |<Soo| 2 versus the scattering energy for the exact calcu- 
lation (solid line) and the adiabatic model (dashed line). 



2M 



+ e w = E 



(21) 



where is the m-th eigenvalue of the Hamiltonian in region % and K$ its corresponding 
centre of mass momentum. 

Similarly, for region i + 1 the wave function will be written as 



jv(*+ 1 )-i 

\^)(E)) = J2 {A { t l) e' lK - +1)R - B ^e +tK - +1>R ) \N^ n). (22) 



n=0 



At R = Ri the part of the wave function associated with the localised state s = i is 
reflected by the wall. Therefore, continuity of the wave function requires that 



(SCLS;Ni\¥^(E)) 



Ri 



0. 



(23) 



By contrast, the rest of localised states are unaffected by this barrier, so we may require 
continuity of the wave function and its derivative for the components s > i: 



(SCLS;Ns\V®(E)) _ = (SCLS; N s\¥ t+1) (E)) 

Fl — Flo 



(24) 
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FIG. 5: S-matrix coefficients for the ground and first excited states. The solid line is the exact 
calculation. The dot-dashed line corresponds to the USA model, with the prescription (|T7|) for the 
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d 

dR 



(SCLS;Ns\ty {i \E)) 



R—Ri 



d 

dR 



(SCLS;Ns\ty {i+1) {E)) 



(25) 



R—Ri 



Equations (l23|) to (l25|) give rise to a set of 2 (~ — % \ -\- 1 conditions for the scattering 
coefficients A$ and B$ [n = 0,...,N^ — 1). When applied to all regions a total of 
(jfj equations with -0 coefficients is obtained. From these, the y coefficients 

are known, as they are determined by the initial conditions. Therefore, there remain 



N ( N 



N 
2 



N 



coefficients to be determined, that coincides with the number of 



2 v 2 1 / 2 V 2 

equations. Thus, the resolution of the problem reduces to the calculation of the inverse of a 
(t) x (t) ma t rix - This can be a slow computational task when iV takes large values. This 
drawback has led us to adopt an alternative method that avoids this sort of calculations, 
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speeding significantly the computational time. For its application, it is convenient to define 
a generalised "S-matrix" for region i as: 



n' \ 



K K 
j\ n , 







^n'n^n' ■ (26) 



K, 



For region i — 1 this definition gives just the usual S-matrix. The factor J K$ / Kn has 
been introduced to ensure unitarity. 

Substituting this expression into equations (123|) to (125]) we get a set of 2 (y — + 1 
equations relating the S-matrix elements of and S^ 1 **. Then, instead of solving the 
system of equations as a whole, we perform an iterative calculation in which the S-matrix 
for each region is determined in terms of the S-matrix for the neighbouring region. The 
starting point is the region i = y, for which the S-matrix is known. In this region there is 
only one state left, which is completely reflected at Rn/2- The S-matrix in this case is just 
the phase factor: = exp ( — 2iK[ 2) R N/2 ). From this, we can determine the S-matrix 

for region y — 1, and so on. Finally, we obtain the S-matrix, in terms of S^ 2 \ Taking 
into account f )26|) and the fact that Affl are given by the initial conditions, it is possible to 
determine the coefficients One can derive also an iterative procedure to determine the 
coefficients A$ for all regions in terms of the corresponding S-matrices. 



C. Discussion of the exact results 



In order to compare both treatments we have plotted in figure [2] the S-matrix coefficients 
<SoO) Sq2 and Sqa, as a function of the scattering energy. The energy scale is in units of Uu. 
A basis space with 15 even parity HO states was used for the calculation. The figure shows 
a good agreement between both treatments, specially at low scattering energies. As the 
scattering energy increases, the effect of the truncation of the basis becomes more impor- 
tant and both calculations differ slightly. Notice that this discrepancy appears to be more 
evident for the phase of the S-matrix elements. Nevertheless, this difference between both 
calculations is reduced as the number of basis states is increased. 

We notice that, as expected, the elastic coefficient is identically one for scattering energies 
below 7:7kij. This corresponds to the energy of the first excited state that is suitable to be 
populated. Accordingly, the inelastic coefficients, Sq2, S04, . . . are identically zero below this 
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threshold. For scattering energies higher than ^hu, the modulus of the elastic coefficient 
is less than one and the rest of the coefficients are non zero. Note that the modulus of 
Sqo presents a minimum at E = 5 hu, indicating a maximal loss of flux from the elastic 
channel to other channels at this energy. For higher scattering energies the modulus of 
elastic scattering coefficient tends gradually to unity. 

We have also compared the wave functions in both approaches. In our problem the wave 
function is a complex quantity depending on two variables, r and R. For the exact direct 
calculation the wave function is given by (jSj), where the scattering coefficients are calculated 
by imposing the boundary conditions ( ITTj) . It must be noticed, however, that this expression 
only holds for n, T2 > or, equivalently, for — 2R < r < 2R. Outside this range, the total 
wave function must be identically zero. The projection of the wave function on a state \Nn) 
is calculated as 

r 2R 

V n (R) = (Nn\V)= dr0* n (r)q(r,R). (27) 

J-2R 

For large values of R (compared to the spatial extension of the internal wave function) 
the integral can be extended to the interval (— oo, oo) and the projection above is directly 
related to the corresponding S-matrix element: 

/oo I 
dr^r^ir, R) = —5 n0 e~ iK ° R - -—S 0n e iK « R ; (R » a ). (28) 
-OO y/V0 y/V n 

In the analysis based on CLS, the total wave function is given in terms of a piecewise 
function of R, according to f[2T)j) . For a certain region i, this wave function is written as a 
superposition of the eigenstates of the Hamiltonian in this region. These eigenstates travel 
freely between the consecutive barriers Ri and Ri+±. Thus, the projection on a state \Nn) 
for a value of R belonging to a region i reads 

atW— l , v 
*®(R) = (iVn|^ W > = E U$e- iK ™ R - B^e +iK ™ R ) (N n\N {i) m). (29) 

m=0 ^ ' 

We notice that, except for the incoming region (i = 1), the states \Nn) and \N^ m) 
are eigenstates belonging to different Hamiltonians and so orthogonality can not be directly 
applied to them. However, they are both given in terms of the SCLS and so the calculation 
of their overlap is straightforward. 
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The results are shown in figure [31 where we present the modulus of the projection of 
the total wave function on the internal states. The solid line corresponds to the treatment 
of section III Al and the dashed line to the CLS method. These wave functions have been 
obtained for a scattering energy E = 5 hu. The R scale is in units of a^. As shown in this 
figure both treatments give almost identical results, the small differences being attributed 
to the truncation of the infinite basis. 

A similar agreement between both approaches is found at other energies, provided that 
a sufficient number of states is included in each case. 

From this picture it becomes apparent the role played by each internal state. The state 
n=0 corresponds to the initial state. Then, the asymptotic wave function contains both 
incoming and outgoing contributions of this state. This is reflected in the characteristic 
interference pattern of the curve | (N 0|\&)|. At this scattering energy, the states n=2 and n=4 
are open channels, and so they give a non vanishing contribution to the asymptotic scattered 
wave. In fact, the curves tend to a constant value for large distances which, according 
( I28jl . is proportional to the corresponding coefficient So n . By contrast, the state n=6 is 
closed. Its exponential decaying tail reflects the fact that this state does not contribute 
directly to the asymptotic wave function, but it does give a non negligible contribution to 
the total wave function in the vicinity of the wall. As mentioned before, the inclusion of these 
states is essential in order to reproduce accurately the boundary conditions and to achieve 
convergence for the S-matrix elements. Then, they indirectly affect the wave function in all 
the space. 

The approach based on CLS presents some advantages compared to the treatment de- 
scribed in section Hi Al It allows to evaluate all the matrix elements S nm , n,m = 0, . . . , N — 1 
in a single calculation, i.e., the initial state does not need to be specified. By contrast, in 
the previous approach a new calculation is required for each initial state. 

Moreover, the method based on the CLS preserves the general properties of the S-matrix 
for any value of N. In particular, conditions ( J23l) and ( J25l) ensure that the total flux is 
conserved at each barrier. As a relevant consequence, the resulting S-matrix fulfils unitarity, 
regardless of the number of initial states chosen for the basis. 
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III. THE ADIABATIC MODEL 



We derive in this section an expression for the S-matrix in the adiabatic approximation 
[lo| . The standpoint of the approach is that a fundamental distinction is made between the 
two relevant coordinates of our problem, namely, the centre of mass coordinate, R, and the 
internal variable, r. The former is identified as a high-energy (fast) variable and the latter 
as a low-energy (slow) variable. 

In this dynamical regime, it is expected that the excitation energies e n associated with 
those excited states which are significantly populated, are such that e n <C E, where E is 
the incident energy of the projectile. Under this assumption it seems reasonable to replace 
the internal Hamiltonian by a representative constant. By choosing this constant as 6q, the 
ground state energy, it is also guaranteed that the solution of the resulting approximate 
three-body equation satisfies the correct incident wave boundary condition. 

Then, applying this approximation to the Hamiltonian Q the Schrodinger equation 
reads: 



+ v l (R + -r) + v 2 (R - -r) + e - E 



q ad (R,r) = 0. (30) 



2MdR 2 K 2 ' v 2 
Note that this approximate Schrodinger equation is independent of the relative momen- 
tum between the fragments. Then, its conjugate coordinate, r, is a constant of motion, 
remaining frozen during the collision. Thus, eq. (1301) has to be solved for all values of a fixed 
separation r. 

In the case of a rigid wall potential, V\ and t> 2 are zero for positive values of r\ and r 2 , 
respectively. Therefore, in this case the solution of eq. (130]) is given by the plane wave, 
exp(— iKqR), multiplied by an arbitrary function of r. The most general solution verifying 
the boundary incident condition at infinity can be written as 

^ ad (r, R) = ct) Q (r)e~ lKoR - S(r)(t) G (r)e iKoR . (31) 

where <po(r) is the ground state wave function and S(r) is a function determined by imposing 
the boundary condition at the wall. This requires that the wave function vanishes at R = 
lrl/2. Then 



S(r) 



-iKo\r\ 
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(32) 



The scattering coefficients defined by eq. (|Sj) can be obtained by projecting the wave 
function fl3Tj) onto the basis states {0 n (?")}• This gives rise to the following simple expression: 



We notice that in this expression no distinction is made between open and closed channels. 
Actually, the approximation treats the full excitation spectrum of the internal Hamiltonian 
as being degenerate in energy with the ground state. As a consequence, unitarity of the 
adiabatic S-matrix is only achieved when summing over the infinite set of eigenstates. 

In what follows we show that the adiabatic approximation is not adequate for the problem 
treated in this work. In Fig. H]the quantity 1 — |Soo| 2 is plotted versus the collision energy for 
the exact (solid line) and adiabatic (dashed line) calculations. This quantity can be inter- 
preted as an excitation probability. The adiabatic prediction completely disagrees with the 
exact calculation, indicating that the assumptions involved in the adiabatic approximation 
are not adequate for this problem. We attribute this failure to the fact that, as revealed by 
the exact calculation, many internal states participate in this process. 

This is probably due to the peculiarities of an infinite zero-range interaction. In this case, 
the momentum transferred to each particle by the wall is twice the incident momentum and 
so, when increasing the scattering energy, the expected excitation energy increases. During 
the collision time, i.e., while one particle has collided with the wall, but the other still has 
not, the internal motion of the projectile is strongly excited. Then, the assumption of the 
adiabatic model, i.e., to consider the whole spectrum to be degenerated with the ground 
state, does not work properly in this case. After the collision, that is, when both particles 
have collided with the wall, the centre of mass momentum is reversed and the final excitation 
energy is small. Thus, this model represents a case in which, although the final excitation 
of the projectile is small, the adiabatic approximation is inadequate, because during the 
collision the internal motion is strongly excited. 

From this discussion we conclude that one should be very careful in applying the adiabatic 
approximation when dealing with strong, short-range interactions. In the next section we 
develop a new method to treat this type of situations. 




(33) 
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IV. UNCORRELATED SCATTERING APPROXIMATION (USA) 



The main goal of this section is to derive an approximated expression for the S-matrix 
corresponding to the one dimensional scattering of a two particle system in terms of the 
constituent-target scattering amplitudes. The results presented here are not restricted to 
the case of an infinite potential. Thus, we start with a general derivation of the model ant 
later we particularize the results to the problem of a wall potential in order to compare with 
the exact solution. 

There are two opposite effects acting on a projectile in the process of the collision. The 
first one is the binding potential v(r) that tends to keep the system bound. The second one 
is the interaction with the target which, apart from governing the motion of the projectile 
centre of mass, is the responsible for exciting or breaking the system. The relative importance 
of these two effects depends importantly on the separation between the projectile and the 
target. In particular, for sufficient large distances between them the dominant interaction 
is clearly the mutual interaction between the constituents. Thus, it seems reasonable to 
approximate the projectile-target potential by an average (folding) potential at sufficiently 
large distances. By contrast, when the bound projectile is close enough to the target the 
dynamical evolution of the projectile is mainly governed by the target interaction. In this 
case it is reasonable to neglect the correlations between the fragments. 

Let us introduce a characteristic centre of mass distance Rq separating these two regions. 
For distances R ^> Rq, referred as the "asymptotic region" and denoted by the index I, we 
adopt an approximate Hamiltonian in which the interaction with the target is neglected 

H^H r + H R = Hi, (34) 
where H r is the internal Hamiltonian and 

Hr = H + V F (R) (35) 

with Vp(R) representing the folding potential between the projectile and target. 

Asymptotically, the eigenstates of Hi are just the product of the eigenstates of the internal 
Hamiltonian H r times a plane wave in R, subject to the restriction As mentioned 

before, there is only contribution to the asymptotic wave function coming from the open 
channels. Denoting by N the number of these states, we restrict the basis space to the set 
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{\Nn); n = 0,1, . . . , N — 1}. Then, the total wave function corresponding to an incoming 
wave in an internal state n will be written in the asymptotic region as 



1 N-l c 

\ty I>n (E)) -^—e- iKnR \Nn) - £ ^ e iK ™ R \Nm). (36) 



For distances R <C R , in what we call "interaction region" (denoted by II), the Hamilto- 
nian (JTJ is approximated by 

H nh 1 + h 2 +v = Hn, (37) 

where v is a constant that substitutes t>(r). 

The eigenfunctions of Hjj can be expanded in terms of the product of eigenf unctions of 
the Hamiltonians h\ and h 2 . An eigenstate of hi corresponds to the distorted wave for the 
scattering of a particle under the potential Vi(ri): 

hi\ X iih)) =E i \ X i{k i ))- (i = l,2), (38) 

where fcj is the asymptotic incident momentum for the constituent i, and E{ = (hki) 2 /2mj. 
Asymptotically this distorted wave behaves as 

|Xi(A*)} ->• l*i> - Si(h)\ - h); (i = 1, 2), (39) 

where (rj|/cj) = exp(— iA^r,)/\/27r is a plane wave with momentum hi and Si(ki) is the 
constituent-target S-matrix for the scattering energy E^. Note that, in the case of the 
infinite wall, this expression is valid for all values of > 0. 

Thus, the eigenstates of the Hamiltonian Hu, corresponding to the asymptotic momenta 
k\ and k 2 , can be expressed as products of the form 

\il>(ki,k2))<x\Xi(ki))\X2(k2)). (40) 

where 

H hi h kr, 

H - + v = E. (41) 

2mi 2m 2 

According to f l39l) . the asymptotic expansion of the eigenfunction (l40l) contains an in- 
coming wave, \ki)\k 2 ), an scattered wave, Si(ki)S 2 (k 2 )\ — — k 2 ), and two cross terms 
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mixing ingoing and outgoing contributions, namely Si(k\) \ — ki)\k2) and S 2 (k 2 ) \ki)\ — k 2 ). 
In practise, these last two terms do not contribute to the incoming nor the scattered wave 
function at large distances. This can be verified by considering an incoming wave packet 
in ki and k 2 . The superposition of states of the form \ki) \ — k 2 ) or | — ki)\k 2 ) contains an 
incoming part which vanishes for t — > +00, and an outgoing part which cancels for t — > — 00. 
Then, these cross terms can be omitted as far as the asymptotic behaviour concerns. As we 
will show later, the use of MLS allows to demonstrate that these terms do not contribute 
to the wave function at large distances in a time-independent formalism. The form of the 
scattered wave, Si^i)^^) | — k\)\ — k 2 ), indicates that the S-matrix for an incoming wave 
with definite values of ki and k 2 , denoted S(ki, k 2 ), appears to be the product of the individ- 
ual S-matrices for the constituents, i.e., S(ki,k 2 ) = Si(ki)S(k 2 ). This S-matrix is unitary, 
provided the individual S-matrices are unitary, i.e., \Si(ki)\ = 1. This condition is satisfied 
for the infinite potential, but it also holds for any other situation for which the transmission 
coefficient is zero. 

Therefore, the scattering problem for the Hamiltonian Hjj corresponding to a situation 
characterised by an incoming wave with definite values of the energies of the constituents can 
be easily solved. However, our physical initial state is not characterised by the individual 
energies of the two particles, but by a certain internal state of the projectile and the energy 
of the collision. The general solution in the interaction region for a total energy E will 
be a certain superposition of eigenstates f j40l) . verifying ( 14TT) and the adequate asymptotic 
boundary conditions. These boundary conditions require that the wave function in region 
II matches smoothly with the asymptotic wave function of eq. ( 1361) . One possible way to 
proceed might be to expand the total wave function in each region in terms of the eigen- 
states of the approximated Hamiltonian for that region. The coefficients of the expansion 
are determined by imposing the continuity of the wave function and its derivative at the 
matching radius Rq. Apart from the complexity of the calculation, this method has the 
problem that incoming waves coming from the asymptotic region do not match exactly with 
incoming waves of the interaction region, due to the discontinuity of the Hamiltonian at Rq. 
As a result, part of the incoming flux is reflected at Rq and spurious outgoing waves are 
generated. To overcome this difficulty we proceed on a different way. In order to avoid the 
unphysical reflections we relax the meaning of the matching radius Rq. We assume that the 
Hamiltonian Hj is smoothly transformed into Hjj in a finite transition region around R . 
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Although we do not make an explicit description of this transition region in our model, we 
include its effect by requiring that the wave function passes from one region to the other 
without loss of flux. 

Let us consider the incoming part of the general solution (1561) : 

|*£?(JB)> -^e-^ R \Nn). (42) 



In order to match this wave function with the inner wave solution it is convenient to 
express the internal states in terms of a basis of Momentum Localised States (MLS). These 
are obtained from the original basis by diagonalizing the momentum operator in the set of 
internal states. Due to the truncation of the original basis the MLS do not have a definite 
value of the internal momentum q but, provided the number of states iV is large, their 
momentum distribution is highly localised around a certain value. The two basis of states 
are connected by an orthogonal transformation (see appendix) and so, the incoming state 
can be rewritten as 

1 N 

\^ m J(E)) -> —— e - iK ^ R ^ (MLS; Ns\Nn) \MLS; Ns). (43) 



'v r , 

where \MLS; Ns) denotes a MLS and (MLS; Ns\Nn) are the transformation coefficients. 
The function (q\MLS; Ns) = (p s {q) has the property of being highly localised around q = q s , 
the s-th zero of the eigenfunction (q\N N). 

As noted before, the eigenstates of the approximated Hamiltonian in the interaction region 
are characterized by the energies of the two particles. These energies are directly related to 
their incident momenta. Equivalently, they can be characterised by the asymptotic values 
of the internal momentum q and the centre of mass momentum K. In our approach, the 
internal eigenstates will be approximated in region II by the discrete basis of MLS. The 
incoming wave function for the interaction region is then expressed at large distances as 

N 

l*/S(-E)> -> E A^e- iKsR \MLS; Ns) (44) 

s=l 

with the centre of mass momentum K s defined by the relation 

^r + w +s = £ - (45) 
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The coefficients are determined in order the waves (1361) and (1441) match smoothly. 
In particular, we impose the incoming flux to be conserved in the transition. This can be 
achieved by taking 



4 n) = -^=e' l5 " s (MLS; Ns\Nn) , (46) 

with v s = hK s /M and 5 ns = 'J n — Xs- The phases 7 n and Xs must be real numbers to preserve 
the incoming flux. 

In order to evaluate the phases 7„ and Xs, we impose the two incoming solutions to have 
the same phase at Rq. This is achieved by taking 7„ = K h Rq and Xs — K s Ro, which leads 
to 

5 ns = (K n - K S )R . (47) 

Once the coefficients are known, the total wave function in the interaction region is 
completely determined. It can be expressed in terms of the distorted waves for each one of 
the constituents. To this end, we rewrite expression (1H1) as 



l&j&m E4 n) e-* K / dq\q)<p a (q) = 

8=1 J -°° 

= J2 A ^ / dqp s (q)\K(q))\k s 2 (Q)), (48) 

8=1 J -°° 

where we have introduced the momenta k{(q) = ^jjK s + q and k^(q) = — q. The 

scattering wave function corresponding to an incoming plane wave \kf(q)) is given by the 
distorted wave \x(k?(q)). Then, the total wave function in region II, including both the 
incoming and scattered wave reads 

\*n,n(E)) = E4 n) / dqUq)\Xi{k{{q)))\X2{k s 2 {q)))- (49) 

8=1 J -°° 

Taking into account the asymptotic behaviour of the distorted waves, (139]) . this wave 
function can be written beyond the range of the potentials as 



9n, n (E!)) -> E4 n) / dqf.{q){\kl(q))\kl(q)) 

8=1 J -°° 
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- S 1 (k s 1 (q)\ - kl(q))\k s 2 (q)) - S 2 (k°(q) |A£( g )> | - k&q)) 

+ ^(^(^(g)! - k{(q))\ - fc-(g))} (50) 

Assuming that the individual S-matrices, Si and 5*2, are smooth functions of the energy 
in the region where the integrand takes significant values, they can be evaluated at k{ and 
k 2 , respectively. Also, it is convenient to express the products of planes waves in terms of 
the relative and centre of mass momenta. Considering the case of equal masses we can use 
the expressions k{{q)r\ + k 2 (q)r 2 = K S R + qr and k{{q)r\ — k 2 (q)r 2 = \K s r + 2qR. Then 
one can perform explicitly the integration with respect to q, to obtain 

N 

*n, n (r,R) -> E4 n) {^V;M + 5iW)5 a fe s )e iif, V s (r) 

s=l 

- S 1 (k s 1 )e i ^ r ( p s (2R) - S 2 (k s 2 )e~^ r V l(2R)} 

(51) 

where <£ s (r) is the Fourier transform of tp s (q). This wave function can be interpreted as 
follows. The incoming wave is decomposed as products of MLS, <£> s (r), describing the internal 
evolution, times an incoming plane wave describing the centre of mass motion, e~ lKsR . This 
wave scatters by the target giving rise to three terms. The second term in (15TT) is just the 
conjugate of the incident wave, times the product of the S-matrices of the constituents. 
The remaining two terms comprise the product of the function (p s (2R), or its conjugate, 
times a plane wave in the variable r. As it can be easily verified, the function ip s (2R) 
vanishes for large values of R and so, these two terms do not contribute to the asymptotic 
wave function. Actually, these terms containing only one of the involved S-matrices can be 
physically regarded as the situation in which only one of the particle has scattered and the 
other has not yet. This is consistent with the fact that they both cancel at large distances. 
However, we remark that these vanishing terms are essential in order to reproduce the 
wave function at small distances. In this sense, it is also interesting to note that the wave 
function (15TT) retains components associated to closed channels, even when they are explicitly 
omitted in the asymptotic region, according ( 1361) . As noted in section II, in the case of the 
wall potential the inclusion of these states is essential in order to reproduce the boundary 
conditions. For this particular problem, expression (l5Tj) is valid for all the interaction region 
and, as can be easily verified, it identically fulfils the boundary conditions, vanishing for 
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r = ±2R. 

Thus, we can conclude that an incoming MLS scatters with the product of the S-matrices 
of the fragments. This result is consistent with our previous discussion in which, using wave 
packet arguments, we concluded that the S-matrix is given by the product of the S-matrices 
of the fragments in the basis characterised by the momenta of the two particles. We now 
see that this property also holds for the MLS basis which, in a sense, can be described as a 
wave packet of plane waves in terms of q, centred around q s . 

Therefore, writing explicitly the value of the coefficients the scattered wave in the 
interaction region behaves at large distances as 



N i 

\^n a n\ E )) ->• E ~7^e~ lSns (MLS; Ns\Nn) e lK ° R S 1 {k s 1 )S 2 {k s 2 )\MLS; Ns). (52) 



N 

By writing the MLS appearing in this expression in terms of the original basis, and 
imposing the conservation of flux, one can easily obtain the scattered wave in the asymptotic 
region 

p iK m R N 
rn \/ v m s= \ 

x (MLS;Ns\Nn}S 1 (k s 1 )S 2 (k s 2 )}\Nm). (53) 

The expression between brackets provides the S-matrix element connecting an initial state 
n with a final state m 

N 

S nm = -Y, e- l{5 ™ +5 ™ ] (Nm\MLS; Ns) S 1 (k s 1 )S 2 (k s 2 ) (MLS; Ns\Nn) . (54) 

s=l 

As pointed before, the peculiarities of the problem treated in this work imply that, 
starting with the system in its ground state, only positive parity states are suitable to 
be populated. As far as expression ( 154"]) concerns, this means that the matrix elements So n , 
with n odd, are identically zero. 

Under these considerations it is possible to exclude those inhibited states. Then, one is 
left with a truncated basis of the form \Nn), where iV is now even. Following the arguments 
above and using the results of the appendix one finally gets the following expression for the 
S-matrix in the case of equal masses 
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S nm = -^ e ~ i(5 '" +5ms) (Nm\SMLS;N s) (SMLS; N s\N n) S l (k s l )S 2 (k s 2 ), (55) 



with n and m even and the sum in s restricted to the positive zeros of "Hjv(x). The coefficients 
(SMLS; N s\N n) are equal to (MLS; N s\N n) up to a factor \/2 (see appendix). 

We remark the simplicity of the expression (1541) as compared to the exact solution of the 
problem. Part of the payoff for this simplicity is the existence of two undetermined param- 
eters, v and Rq. The optimal value of Ro is nevertheless constrained by two considerations. 
On one side, it can not be too large, as at large distances the interaction between the two 
fragments would dominate over the interaction due to the target, and so ignoring the corre- 
lations between the particles for values just below Rq would not be a good approximation. 
On other side, the value of Rq should be large enough to allow us to ignore the interaction 
with the target for R above Rq. Thus, the matching radius must be about the size of the 
system. Moreover, v must be of the order of the expectation value of v (r) on the ground 
state. 

V. COMPARISON OF USA AND EXACT CALCULATIONS 

In this section we analyse the reliability of USA by comparing its predictions with the 
exact results. In Fig. [5] the S-matrix coefficients for the USA model (dot-dashed line) with 
the phases 5 ns given by (14 7p are compared with the exact calculation (solid line) and the 
adiabatic approach (dashed line). The calculations have been performed using the matching 
radius R = 0.6ao and an average potential v = 0. These values were determined by fitting 
the elastic S-matrix at high energies, where the model is expected to be more accurate. We 
notice that for these energy independent parameters a good description of the elastic and 
inelastic coefficients is achieved for energies above ~ 6hu. By contrast, the adiabatic model 
does not seem to give a good description of the S-matrix at any scattering energy. 

A quantity closely related to the S-matrix coefficients is the average final excitation en- 
ergy. It gives an idea on the degree of excitation of the final system. It has been defined 
as 



s 




(56) 



n=0 
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FIG. 6: Average excitation energy. The solid line corresponds to the exact calculation. The dot- 
dashed line corresponds to the USA with the prescription 5 ns = (K n — K s )Rq. 

where the sum must be extended only to the set of open channels. 

In Fig. [6] we show the energy dependence of the average excitation energy, (e). As in 
the previous case, the USA model (dot-dashed line) agrees well with the exact calculation 
(solid line) for energies above 6 hu. On the contrary, the USA model, as expected, does not 
describe properly the low energy regime. In particular, a spurious discontinuity is observed 
for the threshold at E — ~ Tluj. 

In Fig. [7] we compare the wave function given by the USA approach with the exact 
calculation at E — 10 hu. The curve for the exact calculation has been obtained using 
expression ( 1271) . In the USA model, represented by the dot-dashed line, one has to distinguish 
the asymptotic and the interaction regions. For the asymptotic wave function, the interaction 
with the wall is neglected (fj(rj) = 0) and so the range of values of r is unrestricted. Then, 
the projection on a state n is simply given by 



asym 



(R) -> S, 



of) 



-iK R 



0^ e iK n R 



(57) 

By contrast, in the interacting region, explicit account is taken for the wall and an 
expression similar to ( l27j) should be used instead. Thus, the total wave function in II ( I5TT) is 
projected on the different eigenstates, taking into account that the integration is restricted 



27 




1 1.5 

R (a units) 



1 1.5 

R (a units) 



FIG. 7: Projection of the total wave function on the internal eigenstates for a scattering energy 
E = Whoj. The solid line is the exact calculation, the dot-dashed line corresponds to the formula 
(f54j) . with the prescription 6 ns = (K n — K s )Rq. 



to the interval |r| < 2R. 

At this collision energy, Fig. Vindicates that the scattering coefficients are well reproduced 
by ( l54l) with the phases S ns = (K n — K s )Rq. Then, we have adopted this prescription (with 
the same matching radius) to describe the approximated wave function at this energy. As 
can be seen in Fig. [TJ the agreement with the exact calculation is quite good for all the 
states. 

As a general rule, our calculations show a better agreement for the elastic and first excited 
states, and it tends to be worse for excited states of increasing energies. This is expected 
because within the USA model we restrict the basis to the set of open channels and so, the 
effect of this truncation becomes more evident as we explore excited states close to the cut 
off. 

We have explored in more detail the limits of the USA model as well as the validity of 
the prescription ( l47l) for the phases. In order to do that we start from the exact expression 
of the S-matrix for a fixed collision energy. Making use of ( 15^1) . we fit the phases 7„ and \s 
appearing in 5 ns in order to reproduce the exact calculation. We find that the S-matrices 
can be accurately fitted for all the scattering energies (even for very small values). Although 
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these phases, 7 n and Xs, can not be written exactly as K n Ro and K s Rq, respectively, one can 
define effective energy and channel dependent radii so that R n = ^ n /K n and R s = Xs/K s . 
In Fig. [8] we plot the values of the quantities R n versus n and R s = Xs versus the absolute 
value of q s , i.e., the internal momentum at which the state (q\MLS; Ns) is peaked. The 
selected collision energies are E = 9 hu and E = 10 hu. In both cases we deal with a total 
of five HO states, namely, the ground state and the first four excited states with even parity. 
Using the prescription (HT|) both quantities are just Rq for all values of n or s. Note that the 
value of R s is rather constant and very close to the matching radius used in our calculations, 
i.e., Rq = 0.6ao- This constant value has been also plotted in the figure for reference. The 
values of R n are also very close to R = Ro for the lower values of n, but they tend to deviate 
from our prescription for values of n close to the threshold. It is remarkable that the values 
of R n and R s are mostly independent on the scattering energy and of the individual state 
considered. This indicates that for scattering energies large compared to Hoj, the USA works 
very well, and the matching radius can be taken as a constant, related to the size of the 
system, and independent on the energy or the internal state. For lower energies, the USA 
may still be used, but in this regime the radius Rq depends on the energy and the state 
considered. 



VI. CONCLUSIONS 



We have formulated a one dimensional problem consisting on the scattering of two par- 
ticles, interacting with a HO potential, that collide with an infinite potential. This problem 
contains some of the dynamical features of the scattering of composite systems in atomic, 
molecular and nuclear physics, which interact with a target through short range interactions. 

We have obtained the exact solution of the problem using two different procedures. The 
first one consists of imposing the adequate boundary conditions on the scattering wave 
function. The second procedure deals with a basis of configuration localized states (CLS), 
which are wave functions with a strong spatial localisation. Both procedures converge, 
provided that a sufficiently large basis of states is used. The main characteristic of the exact 
solution is that, for large scattering energy, the elastic scattering dominates. In order to 
achieve convergence, the inclusion of closed channels, i.e., states with internal energy larger 
than the scattering energy, is required. Although these states do not contribute to the S- 
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FIG. 8: Behaviour of the ratios R n = ^ n /K n (upper figure) and R s = Xs/K s (lower figure) versus 
the value of n and \q s \, respectively, for the scattering energies E = 9hco (stars) and E = 10?iu; 
(circles). The phases 7 n and \s have been calculated by fitting the S-matrix elements predicted by 
( [54]) to the exact calculation. 



matrix they must be taken into account in the calculations to obtain accurate results. This 
fact indicates that for short-range interactions one should be very careful when truncating 
the basis of states used in continuum discretized calculations. 

We have compared our exact results with the adiabatic approach, that considers the 
relative coordinate frozen during the scattering process. The results disagree completely. 
This indicates that the adiabatic approximation could be inaccurate when the interactions 
of the fragments with the target have a very short range. These short range interactions 
could couple to highly excited internal states for which the adiabatic approximation is not 
valid. However, it should be reminded that our calculations make use of a sharp infinite 
wall. Thus, they will be relevant when the range of the interaction is short compared not 
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only with the size of the projectile but also with the wavelength describing the motion of 
the fragments with respect to the target. 

We have developed a model which describes the scattering of a composite object in terms 
of the scattering wave functions and the S-matrices of the fragments. The model, that we 
call Uncorrelated Scattering Approximation (USA), neglects the correlations between the 
fragments during the scattering. The application of the USA to our model problem gives 
an expression for the S-matrix of the composite system in terms of the product of the S- 
matrices of the fragments. Similarly, the scattering wave function is given as a combination 
of the product of regular wave functions of the fragments. The particular superposition is 
determined by application of the asymptotic boundary conditions. The application of the 
USA to our model problem relies on the use of two parameters. The most important one is 
the distance Rq at which the asymptotic and uncorrelated wave functions are matched. We 
have fixed this value to 0.6a , in terms of the oscillator length, for all our calculations. It 
should be noticed that this parameter has a similar meaning to the matching radius that is 
used in R-matrix theory. The other parameter is the average potential v that replaces the 
interaction between the fragments. We have set this value to v — 0. 

By considering the matching radius to be energy- and state-dependent, the exact S-matrix 
and scattering wave functions are accurately reproduced. Moreover, for high scattering 
energies the elastic and inelastic S-matrices are well reproduced by taking a fixed value of 
the matching radius. 

In general, the agreement with the exact calculation is better for the observables asso- 
ciated with the ground and first excited states and they tend to be worse for states with 
excitation energies close to the total energy. This is attributed to the fact that for these 
states the relative velocity between the fragments is small and then the correlations are 
expected to be more important. 

The one- dimensional model presented in this work can be a useful test case to check 
the validity of different approaches used in the description of the scattering of composite 
systems. The present choice of a sharp infinite wall for the description of the interaction 
with the target, and harmonic oscillator for the interactions between the fragments has the 
advantage that sets fiui as the unique scale for energies and a as the unique scale for lengths. 
In this sense, our results, which are expressed in those units, are valid for any value of the 
mass or harmonic constant. However, the model could be done more realistic, and more 
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complicated, by substituting the infinite wall for an exponential function, and substituting 
the harmonic oscillator by finite potentials. 
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Appendix A 



In this work different sets of basis states have been used, based on harmonic oscillator 
(HO) wave functions. These obey the general form 



0n(r)=A/^exp(-^j-H n Q ; n = 0, 1, . . . (58) 

where = ^Jh/fioj is the oscillator length, M n a normalisation constant and H n the Hermite 
polynomial of order n. In all the calculations we restrict the infinite set of states to a finite 
family of N states, denoted by {\Nn); n = 0, . . . , N — 1}, where (r\N n) = (j) n (r). 

By diagonalizing the position operator in the truncated basis of HO wave functions, a new 
set of N states is obtained which have the property of being localised in configuration space. 
They are called Configuration Localised States (CLS) and are denoted by {\CLS; Ns); s = 
1, . . . , N}. These two families of states are related by means of the orthogonal transformation 



N-1 



\CLS; Ns) = (Nn\CLS; Ns) \Nn) (59) 



n=0 



For HO wave functions, the state (r\CLS; Ns) is localised around r = aox s , where x s is 
the s-th zero of the Hermite polynomial Hn(x). In this case, the transformation coefficients 



are given by (see 22j |) 



(Nn\CLS; Ns) 



l N - n (N - 1)\] 1/2 Uj.r» 



T^-n-i{x s 



(60) 



2N n\ 

In analogy with the CLS, it is possible to define internal states localised in momentum 
space, known as Momentum Localized States (MTS). This is carried out by diagonalizing the 
momentum operator in the truncated HO basis. Thus, starting with the basis of N states, 
this procedure provides a new set of N internal states, each one of them is peaked around a 
certain momentum. As in the case of the CLS there is an orthogonal transformation relating 
both sets of states: 

N-1 

\MLS;Ns) = Y, (Nn\MLS;Ns)\Nn), (61) 

n=0 

where \MLS; Ns) represents a MLS. Then (q\MLS; Ns) is localised around q = x s /ao 
where, due to the formal analogy of the HO wave functions in momentum and configuration 
space, {x s } are again the zeros of the Hermite polynomial Hn- This analogy provides also a 

33 



simple relation between the transformation coefficients (MLS; Ns\Nn) and {CLS; Ns\Nn). 
In this work, we take the internal wave function to be real in configuration space, and so 
they will be affected by a factor (— i) n in momentum space. Moreover, we take the CLS to 
be real functions, so the coefficients (CLS; Ns\Nn) are real numbers. Therefore, 



{MLS; Ns\Nn) = (-l) n (CLS; Ns\Nn). (62) 

Although the formalisms presented in this work do not require the fragments to have 
equal masses, we have performed all the calculations for this particular situation. In this 
case, parity conservation implies that only positive parity states are suitable to be populated 
during the process. In fact, the coefficients So n derived from the exact calculations of sections 
(III Al) and (III Bl) are found to be zero for odd values of n. This is also satisfied by eq. (154"1) . 
as can be easily verified. 

Under these considerations one is allowed to exclude those inhibited states from the 
beginning. This permits to work with the truncated basis of states: {<fi n (r); n = 0, 2, . . . , N— 
2}, with iV even. This requires, however, some care in the evaluation of the localised states 
in the new basis. The formalism of CLS can not be directly applied to this set, as many 
of its properties entails the sum over both even and odd states. The star ting point for the 



construction of the CLS formalism requires a set of functions of the form 



22J 



ip m (x) = (x\jm) =N jrn \ /2 F{y)V m {y), m = 0, 1, . . . , j - 1, (63) 

where y is a certain function of r, J\fj m a normalisation constant, F an arbitrary function of 
y and V m is a polynomial of order m. 

The drawback outlined above can be overcome in the case of HO wave functions by 
writing the Hermite polynomials in terms of generalised Laguerre functions (see, for instance, 
Ref. 123(1): 

fo m (r) = (-l) m m!2 2m Ad exp C^iv), V = {^) • (64) 

Then, taking V m (y) = £ { m 1/2) (y), F(y) = exp(-y/2) and j = f we can identify if) m (x) = 
02m(^) (m = 0, 1, . . . , y — 1). The set of configuration localised states are now calculated for 
the new set of functions. It requires the calculation of the roots of the polynomial C^H^ (v) 
which, attending to (164]) . are just the square of the zeros of 1-Ln{x). The new set of localised 
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states are given by a linear combination of the states n (r) and so they are even functions 
with respect to the variable r. Therefore, we call them Symmetric Configuration Localised 
States, SCLS. The transformation between the set of states (151)1 and the SCLS is then 
expressed as 

^ N 
\SCLS;Ns)= Yl (Nn\SCLS; Ns)\Nn), s = l, (65) 

n=0(even) " 

The state (r\SCLS; Ns) is localised around r = ±aoX s where {x s } are the positive zeros 
of l-L N (x). 

The coefficients (SCLS; N s\Nn) are found to be equal to those appearing in eq. 



up to a factor \/2: 

{SCLS;N s\Nn) = V2{CLS;N s\Nn) (66) 

for n even. 

As an example, in Fig. Othe set of HO wave functions n = 0,2,4,6 are plotted versus 
the adimensional variable x (upper figure). The corresponding Symmetric Localised States 
are also plotted (lower figure) and labelled with the index s (s = 1,2,3,4). Notice that 
each one of these localised states is peaked around two symmetrical points, corresponding 
to symmetrical roots of the Hermite polynomial Tisix). 

In a similar way, it is possible to construct Symmetric Momentum Localised States, which 
are given by means of the transformation 

*Z=? N 
\SMLS;Ns)= ^ (Nn\SMLS; N s) \Nn), s = 1, — . (67) 

n=0(even) ~ 

The state (q\SCLS; Ns) is localised around q = ±x s /a where {x s } are again the positive 
zeros of Hj^(x). In the case of the HO basis, the transformation coefficients are related to 
those in configuration space: 

(SMLS;Ns\Nn) = V2(-i) n/2 (CLS; N s\N n) . (68) 
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